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The perpendicular magnetic anisotropy (PMA) arising at the interface between ferromagnetic 
transition metals and metallic oxides are investigated via first- principles calculations. In this work 
very large values of PMA up to 3 erg/cm 2 at Fe|MgO interfaces are reported in agreement with recent 
experiments. The origin of PMA is attributed to overlap between 0-p z and transition metal d z i 
orbitals hybridized with d xz ^ yz ) orbitals with stronger spin-orbit coupling induced splitting around 
the Fermi level for perpendicular magnetization orientation. Furthemore, it is shown that the PMA 
value weakens in case of over- or underoxidation when oxygen p z and transition metal d z i orbitals 
overlap is strongly affected by disorder, in agreement with experimental observations in magnetic 
tunnel junctions. 

PACS numbers: 75.30.Gw, 75.70.Cn, 72.25.-b, 73.40.Rw 



Spin-orbit interaction (SOI) plays a major role in a 
wide class of physical phenomena both from fundamen- 
tal and applications points of view For instance, it is 
at the heart of basic magnetic phenomena such as mag- 
netocrystalline anisotropy Rashba effect 0, Q and 
magnetization damping. Controlling SOI strength at the 
interface between ferromagnetic (FM) and non-magnetic 
layers represents an outstanding challenge for advance- 
ment of transport and magnetic properties of spintronic 
magnetic devices, such as perpendicular Magnetic Tun- 
nel Junctions [5Ml0j (p-MTJs) and tunneling anisotropic 



magnetoresistive (TAMR) systems [ll|. Recently, elec- 
tric field control of interfacial magnetic anisotropy has 
attracted much attention as well 12, [Uj. Traditionally, 
interfaces between magnetic and heavy non- mag netic 

transition metals such as Co|Pt 0, co|Pd Ha B3, 

Co|Au [13] have been used to obtain perpendicular mag- 
netic anisotropy (PMA). It has been shown that the onset 
of the PMA at these interfaces is related to an increase 
of the orbital momentum of Co [13] due to the strong 
hybridization between the 3d orbitals of the transition 
metal and the bd orbitals of heavy metal [l4| . This hy- 
bridization enhances the energy splitting between the Co 
3d z 2 and 3d x i_ y i orbitals and induces a charge transfer 



between the two layers [18M2CH] . As a result, the com- 
bination between SOI and hybridization- induced charge 
transfer leads to the PMA. Thus, the presence of a heavy 
non-magnetic layer (Pt, Pd, Au, W, Mo) was believed to 
be essential to trigger the PMA. 

However, Monso et al have shown that PMA could 
be observed also at Co(Fe)|MOx interfaces (M=Ta, Mg, 
Al, Ru etc.) [21] in spite of the weak SOI at the inter- 
face. Surprisingly large PMA values up to 1 to 2 erg/cm 2 
have been reported, which are comparable or even larger 
than the PMA observed at Co|Pt or Co|Pd interfaces [221 ]. 
This result is quite general and has been observed in both 






FIG. 1: Schematics of the calculated crystalline structures 
for (a) pure, (b) over-oxydized and (c) under-oxydized ge- 
ometries. Fe, Mg and O are represented by blue, green and 
red balls respectively. 



crystalline (MgO) or amorphous (AlOx) barriers, using 
both natural or plasma oxidation [23|. The PMA could 
be dramatically improved under annealing [To|, an d X- 
ray Photoemission Spectroscopy has demonstrated that 
the PMA could be correlated without ambiguity with the 
presence of oxygen atoms at the interface [23|. In fact, a 
correlation between PMA and oxidation conditions have 
been demonstrated for a wide range of FM|MOx includ- 
ing those based on Co^Fei-^, thus indicating that the 
phenomenon is quite general at interfaces between mag- 
netic transition metals and oxygen terminated oxides. 
These observations led the authors to postulate that, in 
spite of the weak SOI of the elements (Fe, Co, Al, O), 
oxidation condition plays an essential role in the PMA, 
as it does in TMR [25| or interlayer exchange coupling 
(IEC) [26j. Recent experiments reported large PMA val- 
ues of 1.3 erg/cm 2 at CoFeB|MgO structures 0, EH- 



2 



-74.3349 
-74.3350 




30 60 
9 (degree) 



FIG. 2: Angular dependence of the magnetic energy, where 
is the angle between the magnetization direction and the 
normal to the interface plane. 



Furthermore, it has been experimentally demonstrated 
that there is a strong correlation between PMA and TMR 
maximum values obtained at the same optimal oxidation 
and annealing conditions [29| . 

In this letter, we report first-principles investigations 
of the PMA and the effect of interfacial oxidation condi- 
tions on the PMA at Fe|MgO(100) structures. The latter 
can be viewed as a model system for FM|MOx interfaces 
involving bcc electrodes including Co^Fei_ x alloys. In 
agreement with experiments, it is demonstrated despite 
the weak SOI, the bonding between the Fe-3d and O- 
2p orbitals gives rise to even stronger PMA compared to 
that of Co|Pt interfaces. The largest PMA value is ob- 
tained for ideal interfaces while it is reduced for the case 
of over- or underoxidized interfaces. 

For the ab inito calculations, Vienna ab initio simula- 
tion package (VASP) J30| was used with generalized gra- 
dient approximation [31( and projector augmented wave 
potentials [HI . The calculations were performed in three 
steps. First, full structural relaxation in shape and vol- 
ume was performed until the forces are smaller than 0.001 
eV/A for determining the most stable interfacial geome- 
tries. Next, the Kohn-Sham equations were solved with 
no spin-orbit interaction taken into account to find out 
the charge distribution of the system's ground state. Fi- 
nally, the spin-orbit coupling was included and the total 
energy of the system was determined as a function of 
the orientation of the magnetic moments. A 19x19x3 
K-point mesh was used in our calculations with the en- 
ergy cut-off equal to 520 eV. Three structures were con- 
sidered as shown in Fig. [TJ (a) "pure" (O-terminated) 
interface, (b) over-oxidized interface (with O inserted at 
the interfacial magnetic layer), and (c) under-oxidized 
(Mg-terminated) interface. We point out that the situa- 
tion of "pure" interface is the most stable one as observed 
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FIG. 3: Spin-orbit coupling effects on wave function char- 
acter at T point of interfacial Fe d and neighbor oxygen p z 
orbitals for pure Fe|MgO interface shown in Fig. [TJi. Three 
subcolumns in each column show the band levels for out- 
of-plane(left) and in-plane(right) orientation of magnetiza- 
tion as well as for the case with no spin-orbit interaction in- 
cluded(middle), respectively. Numbers are the percentage of 
the orbital character components within Wigner-Seitz spheres 
around interfacial atoms. 



in annealing experiments [24[. The most stable location 
for the oxygen atoms is on top of metal ions due to strong 
overlap between Fe-3d and 0-2p orbitals. Correlatively, 
it is intersting to note that this structural configuration 
also yields the spin filtering phenomenon based on Bloch 
states symmetry leading to high TMR values [HI, liij . 
Furthermore, the strong hybridization significantly mod- 
ifies the band structure giving rise to a high interfacial 
crystal field [H|. 

In Fig. [2j we present the calculated energy per unit 
cell as a function of the angle between magnetiza- 
tion orientation and normal to the plane for pure un- 
relaxed Fe|MgO interface. The dependence is well fitted 
by the conventional uniaxial anisotropy formula Ea = 
K + K 2 sin 2 6>, where K 2 = 0.7 meV/atom (K 2 = 
1.36 erg/cm 2 ). Interestingly, the perpendicular surface 
anisotropy in this case is stronger than that of Co|Pt [22| 
in agreement with recent experiments 27]. The calcu- 
lated anisotropy value is further enhanced for relaxed 
structures as shown in Table HI The PMA for relaxed 
structures weakens in case of interfacial disorder and be- 
comes equal to 2.27 and 0.93 erg/cm 2 for under- and 
overoxidized cases, respectively (see Table I]), indicat- 
ing that the oxidation condition plays a critical role in 
PMA as it does in both TMR ^ and IEC [26]. Fur- 
thermore, the tendency of PMA to decrease with oxy- 
gen excess or deficit along the metal/oxide interface is 
consistent with the recent experimental observations of 
PMA dependence on annealing temperature and oxyda- 
tion conditions Io|, 29| . It was reported that with higher 
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FIG. 4: The same as Fig. [3] for over-oxidized Fe|MgO inter- 
face. 



FIG. 5: The same as Fig. [3] for under-oxidized Fe|MgO inter- 
face. 



annealing temperatures PMA increases due to interfa- 
cial quality improvement (Ioj |. As we stated above, PMA 
reaches a maximum value corresponding to the TMR ra- 
tio maximum indicating that ideal interfaces are crucial 
also for PMA observation [29j. 

In Table U we give also the evolution of interfa- 
cial Fe magnetic moments as a function of distance 
from the interface. One can see that compared to the 
pure case, the moments are enhanced (weakened) for 
over (under) oxidized interfaces. 

Let us now proceed with the explanation of the phys- 
ical origin of the results obtained from first principles 
on the effect of oxidation conditions on PMA. To un- 
derstand the PMA origin at Fe|MgO interfaces, we per- 
formed detailed analysis of the impact of spin-orbit in- 
teraction on electronic band structure with out-of-plane 
(d z 2, d XZl d yz ) and in-plane (d x 2_ y 2, d xy ) Fe-3d and 0-p z 
orbit als character. 

We start from the analysis for pure interfaces repre- 
sented in Fig. [TJl. In Fig. [3j we show bands around the 
Fermi level Ep at T-point with orbital and interfacial 
atoms projected wave function character for out-of-plane 
(left) and in-plane (right) orientation of the magnetiza- 
tion as well as in the absence of spin-orbit coupling (mid- 
dle). When no SOI is included (middle subcolumns), one 



TABLE I: PMA(in unit of erg/cm 2 ) and magnetic moment 
(Hb per Fe atom) for different layers of Fe in Fe|MgO MTJs 
with different oxidation conditions. 
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2.93 


2.27 


0.98 




interfacial 


2.73 


2.14 


3.33 
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sublayer 


2.54 


2.41 


2.70 




bulk 


2.56 


2.55 


2.61 



can clearly see the band level resulting from hybridization 
between Ye-d z 2 and 0-p z orbit als. This is a signature of 
Bloch state with Ai symmetry around Ep for majority 
Fe and MgO which is at the heart of the spin filtering phe- 
nomenon causing enhanced TMR values in MgO based 
MTJs [33|. The double degenerated bands with A5 sym- 
metry related to minority Fe are also present close to the 
Fermi level. When spin-orbit interaction is switched on, 
the degeneracy is lifted and majority Ai and minority 
A5 are mixed up producing bands with both symmetry 
characters. As a result, band levels with d z i, d XZl d yz 
and p z character appear splitted around the Fermi level 
and this splitting is larger and the lowest band deeper for 
out-of-plane magnetization orientation as clearly seen in 
left and right subcolumns in Fig. [3j respectively. Thus, 
the lift of degeneracy of d xz and d yz orbitals combined 
with mixing with Fe-d z 2 and 0-p z orbitals is at origin 
of perpendicular magnetic anisotropy for pure Fe|MgO 
interfaces. This result shows that the out-of-plane com- 
ponents of d xz ^ yz orbitals plays a crucial role for PMA 
similarly to Co|Pd interfaces [35]. 

Next, we proceed with the same analysis for under- and 
over-oxidized Fe|MgO interfaces represented in Fig. [T](b) 
and (c), respectively. As shown in Fig. 2] for the case 
with an additional oxygen located at the Fe|MgO inter- 
face (Fig. [H(b)), spin orbit coupling lifts again the de- 
generacy for states with d xz ^ yz causing stronger splitting 
and deeper level position for out-of-plane orientation of 
magnetization compared to the in-plane one. However, 
these states are not mixed anymore with Fe-d z 2 and O- 
p z orbitals due to local charge redistribution induced by 
additional oxygen atoms [25]. Since d z 2 and p z orbital 
hybridization which is mainly responsible for PMA is not 
splitted, the anisotropy is significantly reduced. 

A different picture occurs in case of underoxidized 
Fe|MgO interface represented in Fig. Die). As shown 
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FIG. 6: Ai Bloch state character at T-point around the Fermi 
level as a function of layer number in pure and over-oxidized 
Fe|MgO interfaces shown in Fig. [T^a) and (b), respectively. 
Ai Bloch state is absent around Fermi level in underoxidized 
case shown in Fig. [TJc). 



in Fig. [5j the Fe-d z 2 and 0-p z components around the 
Fermi- level are now absent. As a result, the degener- 
acy lift induced by spin-orbit interaction for states with 



^xz^yz 



character is now solely responsible for the PMA. 
Since the splitting of these d XZjyz orbitals is still relatively 
strong and higher for out-of-plane magnetization orienta- 
tion compared to the in-plane one, anisotropy values are 
higher compared to overoxidized case but lower compared 
to the ideal Fe|MgO interfaces. Thus, the PMA reaches 
its maximum for ideal interfaces. In order to understand 
the correlation between PMA and TMR, in Fig.[6]we plot 
the wave function character of Ai Bloch state as a func- 
tion of the position across the supercells used for PMA 
calculations. One can clearly see that the Ai decay rate 
is strongly enhanced in the case of overoxidized interface 
compared to the ideal one. There is no Ai band around 
the Fermi level for underoxidized case as demonstrated 
above. This explains why both PMA and TMR reach 
the maximum values in a correlated way as observed ex- 
perimentally [29], this maximum being reached for ideal 
interfaces. 

In conclusion, we presented ab initio studies of per- 
pendicular magnetic anisotropy at Fe|MgO interfaces as 
a function of the oxygen content along the interface. The 
PMA values are highest in the case of pure interfaces in 
agreement with recent experimental studies 27, 2! 
may reach up to 3 erg/cm 2 for relaxed interfaces 



and 
The 



origin of large PMA is ascribed to the spin-orbit induced 
mixing between Ai and As-like 3d and 2p orbitals at the 
interface between the transition metal and the insulator 
combined with the degeneracy lift of out-of-plane 3d or- 
bitals. The PMA amplitude degrades in the case of over- 
or underoxidized interfaces in agreement with recent ex- 
periments [Io|, Hsj. This is due to the impact of split- 



ting(disappearing) of Ai-like hybridized states around 
the Fermi level in presence (absence) of additional oxy- 
gen atom. 
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